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^l2(xrx2)  ~ r]_Cx]_)p2  Cx2^1  + aSl  ^xl^32 


(i/itli  | ot  | < I ) 


where  F.  (x.)  is  the  cdf  or  X.  (j=l,2)  and  S.  (x.)  = 1 - F, (x.)  is  tne 

J » J J J J J 

survival  distribution  function  (sdf)  of  'A-  . 

In  Johnson  and  Xotz  (1976)  vie  discussed  certain  generalizations  of  (1) , 
Motivated  by  sone  remarks  of  Farlie  (I960),  we  suggested  replacing  the 
product  5-l  (ic^ ) S2 (ec0)  (which  would  be  the  joint  sdf  of  X,  and  X2  if 
they  were  independent)  in  (1)  by  the  joint  sdf  of  another  bivariate  FGj 
distribution,  namely 

^12^  (xl»x2)  = (^2)  ^1  + ot^r*^Cx^)F2(x2)  1 (2) 

with  | | £ 1 . 

Tnis  procedure  can  be  iterated  in  various  ways.  In  particular,  we  cau 
replace  in  (2)  by  the  FG:  type  distribution 

F12)(x1’x2)  = F1(x1)F2(^U  +a2Sl(V32(x2)}-  (3) 

(2) 

(Of  course  Fj-ip  might  itself  be  <y‘  the  foxn  of  a -irst  stage  iterate 
FGl  distribution,  like  (2),  but  here  we  will  restrict  ourselves  to  the 
relatively  simple  form  (3).) 


One  purpose  of  introducing  these  iterated  generalizations  of  bivariate 
FGP  distributions  was  to  check  if  one  could  remedy  a defect  of  bivariate 


2 


FCTI  distributions  — namely,  the  relatively  low  correlation  coefficients 
tiiat  can  be  attained  with  such  distributions. 

For  FG!  distributions  with  nomal  marginals  the  correlacicn  coefficient 
-» 

is  air  A , and,  since  |a|  s 1,  this  neons  that  the  absolute  magnitude 
of  tlie  correlation  coefficient  cannot  exceed  tt  ~ (=0.32  approx.) 

In  the  present  note  we  find  that  iterations  of  the  kind  described 
above  cannot  increase  the  maxima  possible  dependence  (as  measured  by  an 
index  suggested  by  Schweizer  and  V/olff  (1976)),  between  the  variables  by 
more  than  about  0.12,  no  natter  how  many  iterations  are  performed,  ’/e 
therefore  investigate  a few  other  nodifications  of  bivariate  FGL!  distribu- 
tions which  night  increase  dependence  more  substantially,  but  vie  find  that 
necessary  restrictions  on  the  values  of  the  ads  are  such  as  to  seriously 
limit  the  possibilities. 


1.  fcdweizer  6.  Wolff's  [iessure  of  dope  no  e nee 

The  measure  of  dependence  introduced  by  Schweizer  § rolff  (1976)  is 
essentially  (for  continuous  joint  distributions) , the  expected  value  of 


12|F12CdrX2)  - F1(X1)F2(X2) 


(4) 


(or;i  equivalently,  of 


i2|s12(x1,x2)  - s1(x.1)s2(r.2)|  ) 


(4’) 


Denoting  this  by  the  symbol  (aw)v  v it  is  clear  that  it  is  unchanged 

'VI 

by  any  monotonic  transformation  is  applied  to  either  (or  botn)  of  the  variables 


XlfX,  . 


In  particular  we  can  suppose  that  such  transformations  have  been 


■5 


applied  to  produce  variables  Xp  X,  which  are  each  uniformly  distributed 


between  0 and  1.  Then 


(crw)v  v 


f1  f1 

12  |F17(xn,x^)  - x,x?|dx,c 

Jn  Jn  LL  1 “ 


If  X,  and  X,  are  mutually  independent,  (ow)v  v = 0.  The  maximum 

'VJ'-2 

possible  value  of  (aw) ,,  Y is  1,  arid  this  is  attained  when,  with  proba- 

a1,a2 

bility  l,  Xj  arid  X2  are  nonotonic  functions  cf  each  other 

Chen  Xj  and  ^ are  each  uniformly  distributed  over  (0,1),  if 

X^  and  7-2  are  each  rronotonic  functions  of  the  other  then  we  rust  have 
either  X^  = X~  or  X,  - 1 - X0  .) 

he  will  use  this  measure  of  dependence  in  our  studies,  rather  than 
the  better- known  mean  square  contingency.  This  is  quite  closely  related 
to  the  Schweizer-  'olff  index, (see  Section  5)—  but  the  comparisons  using 
it  are  not  so  clearcut  as  using  the  latter  measure. 


3 . deea.v-’ence  in  Generalized  fh:  distributions 

If  tie  marginal  distributions  are  each  uniform  over  (3,1)  then  the 


FJSI  cdf  (1)  is 


x:x2(l  + a(l-xi)(l-:c2)} 


rl  (1 

(oio)Y  y ~ 12) dt | x,  (1  -x.  )xhl-x.)  nc,  ix? 

AlAi  Jq  Jr)  x ± *•  ~ 

- 12 |ot|  (j  - \)2  = ||a|.  (0) 

T!iis  is  o(  course  the  value  of  (aw)v  Y for  any  continuous  FfT  • distribution 

u'i2 

v/ith  cdf  (1) . 


f;  y k b I 'Lr*1 


> Am 


l 
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For  the  once  iterated  FG»’  distribution  vdth  unifom  (0,1)  uarginal 
distributions  we  have 


F12^x1^x2^  = xix2^-1  + a(l-xi)(l-x2){l  + oijX^}] 


whence 
(oui)  , „ = 

At  jV, 


rl  rl  - _ 

y = 12  |cm,  (l-x,)x-,(l-x.,)  + aa,x'(i-x,)x^(l-x,)  |dx,dx?  . 

'^'2  J o t Q * X it  L X X XL  X XL 


For  51 ven  numerical  values  of  a,  0^  we  v/ill  get  the  largest  value  for 
(ou))v  by  having  a,,  a,  of  this  sane  sign.  Taking  a>3  and  a,  > 0 

X X 


we  oocam 


C^X,,.  = 12 {a (j  - t)“  + aa1(i  - } * ia  + ^aon  « i(l  + ^a^a.  (3) 

Similarly,  vs  find  for  the  twice  iterated  FQ<i  distribution,  with 

a ^ 0,  a.  £ 0.  a-  £ 0 


= b + i7°ai  + 7dciaict2 


Genera] ly  for  the  r-th  iterated  i'CT ! distribution  with  no  negative 


where  aj  = a 


(aio).-  y =12  X a,a! 
1*2  j=--l  J J 


3-1 

a!  = a T!  a (j*2) 
J i=l  i 


for  j odd 

,n(Xp;  Xp)  }2  for  j even 


where  B(i,j)  is  the  Beta  function. 


Sinc^  Osa^s^i  ...  iaJ«o,  it  follows  that 


and  so 


(aio)v  v s 12a  t a.  < 12a  Y a. 
h*2  j=i  j j=i  J 

(ow)Y  v ^ 12a  l [{B(h,h)}2  + {B(h+ls  h)}2l 
X1A2  h=2 

oo 

= 1 5a  l {3(h,h)  Y" 
h=2 

(since  501+1,  h)  = %B(hsh)  ). 


iov7  l (3(h,h)}2  - 6'2  + SO"2  + 140'2  + G30'2  + ... 

2 

< 6"2  + 30"2  + 149' 2 + 650"2(1  + 4'2  + 4-4  + 


and  so 


< 0.029 


(au>)Y  Y < 15  x 0.022a  = 0.435a 

' j.  2 


however  nan/  iterations  are  used. 

Tins  snows  tiiat  iterations  of  the  kind  described  in  Section  1 cannot 

increase  the  measure  of  dependence  (aw).,  above  0.435. 

'“'‘1 } 

Tiie  reason  for  the  low  limit  on  the  measure  of  dependence  is  in  the 
rapid  decrease  of  B(h/n)  as  h increases.  A natural  way  to  avoid  this 
would  be  to  use  distributions  of  forr, 

F12(xl*x2)  = F1(X!^F2^X2^  '-1  + ct'f5^(x1)}  1{S2(x2)}<^>2] 


with  0 < 4>  - < 1 (j  = 1,2).  Unfortunately  this  is  not  a proper  distribution 
function  for  any  absolutely  continuous  marginals  Fjfxj),  F2(x2).  Ve  have 

fjCxjjfjfctjjri  . , «,(*,))■  ‘b  U - (wp^txp)]. 


Cs Ji  I 


2y  dwosii'c  '-'i  so  that  (1  + ♦1)F(x1)  > i sad  so  that  S^Cx.J  is 
sufficiently  snail  we  tta.-ce  hie  e:c  ression  in  severe  Lractets  no "Stive. 
(Since  a (x,)  < 1 aiti  2 (i^)  > 0,  ire  can  alvrays  find  x-, , x..  so  that 
£-  (x-; ) > G and  £7(x->)  > 0 .) 

— CJ  L* 

3os  ie  otner  ratifications  of  bivariate  rVl  .’istributicns  are  diseased 
in  t ■£,  next  section,  witn  special  reference  to  nossible  decreases  in  deven- 


*. . tone  vtber  Oa.'ijirti  itsd  f-V  ,-istrl1. jfiuQi 

Yne  folio  .'in-  exsirryle  shov.'s  directly  her*?  die  restrictions  on  she  values 
of  the  a’s  militates  against  realization  of  substantial  increases  in 

(VC)j.r  v . 

*ll42 

"ro  consi  der  a faroiiy  of  distributions  obtained  by  rodiryirK;  (1)  vith 
re -lacs:  lent  of 


h^-2^2)  by  Fi(-i)F2(x2){i  + c^Cx^S^)} 
and  C,  (x  )G  ,(m7)  by  S.  (x,)S2(x;)  {1  + a2F,  (x,  )F.  (r'} }. 

7r*uis£6min£  so  tliat  each  nsryinsl  '.istribution  is  uniform  over  (0,1) 
ve  have 

12^\if'tP  = + al (l-^)  Hu  + rt(l~Xj)  (lrx7)  d + a^x,)) 

(0  < 1,  j=l,?.).  (15) 


^X, 


r1  r 1 ? 

" 12  j.  J,  Ka+a^x^d-xOd-;^)  + cca,^ (!-:<,) 

+ aa-nc^x. (l-r.i) (l-x7)  + o^ayy/i (l-/^) J (l-n7) ~ | nr-, con 


iiuii  hiAi*!* WkL _CU(1 


L >j  to 


If  all  a'Js  are  norms  - at  ive  then 


(ow). 


1 1 -1 
^ v = -r( a + a., ) + ipa  (a,  + a-, ) + yr-acxja^  . 


(16) 


If  we  only  needed  the  conditions  |c;|  < 1,  laj  < 1,  |a0|  < 1 then  it 

would  appear  that  we  could  set  values  o£  (ow)-  v as  near  to  - 

11  ' -VA2  - • 

+ g-  + jtj  = 0.346  as  desirea. 

However » the  a’s  must  be  such  tliat 

7 

Txfaij=  1 + (a+di)  (l-2x1)  (1-2x2)  + m1(l-x1)(l-3x1)(l-x2)(l-3x2) 

+ oa2xi(2-3x1)(2-3x2)  + 4aa^a2x,  (1-Xj)  (l-2^)x2(l-x2)  (l-2x2)  > 0 (17) 


for  Xj,  x2  with  0 s xi  s i,  0 s x2  s 1, 
Piloting  = 0,  x2  - 1 gives 


a + a-,  si. 


This  implies  that 


(au))H.,A„  - I + 12°(1  + a2  • a>  + 7?*V-a)«2 

1 L, 

5 T + X2a(2'a)  + 75«(l-a) 


s 0.417  for  Osasl. 


(13) 


fhe  modified.  FG  ••  fai  lily  defined  bp 


r12^xl,x2^  = ^l^x1^2^x2^  ‘ ^ + a C^j)  > S?(x2)Jj 


(19) 


or  equivalently 


\~%s 


F12^xl,x2^  = WW*  + * rax-(F,  (x,).F7(x7))}]  (19)  ’ 


has  a somewhat  larger  value  of  (oo>)  - 

rl  rl 


V'2 


x - i-  ■ 2 l 2J 

We  have 


(ow)v  v = 12  a 
*1*2 


0 JO 

= 12  a x 2 


x.x7{l  - max(XpX2)  }dx^dx2 

i X2 

L L 


= 12  a 
= 0.6a. 


.3 


X2(l-X2)dx2 


(20) 


This  is  considerably  larger  (if  a>0  ) than  the  values  so  far  obtained, 
but  it  still  cannot  exceed  0.6. 


5 . Comparison  with  ilean  Contingency 

The  mean  contingency  (as  defined  by  R6nyi  (1970,  p.  232)  ) for  an 
absolutely  continuous  distribution  with  standard  uniform  marginal  distribu- 
tion is 

? r1  r1  -> 

(21) 


32F 


x2 

< p-r  yr 

A1}A2 


0 


' Cf 12 .Xj) 32dx1dx2  • 1 
0 


12  • 


where  = JxTx~  Fs  t'ie  probability  density  function  of  and 

h • 

Table  1 summarizes  values  of  (ow)Y  Y and  the  corresponding  values 

2 '1^ 
of  <j)v  ,,  , for  the  distributions  discussed  in  tiiis  note. 

*1*2 
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TABLE  i Comparison  of  (ay) Y v and 


u>  , 

Value  of  (au)„  v 
a1A2 

Corresponding  value 

7 

from  ea. 

value 

of  <r 

(6) 

ji“i 

xxx 

12  ^ 
9° 

(3) 

^1+Ictl+I5al)a2 

0) 

3a4zactl+T5ctctla2 

In.l  4 2 2 2 2 4 2 2,  2 

7(1+2^1+23al+2^1a2+2^1a2+Tm-ala2)a 

(16) 

-‘-(cx+a-j^)  +JJ01  (a1+a2)  +y-cta,  a2 

g-{  (a+o.1)2+^a(a+a1)  (a1+a2)+^a2(a^+at) 

+t^(-a+al^aaIa2+575a  a^G12^0lI+a2+^^ 

4 2 2 2, 

ala2} 


note  that  for  (6) 


and  for  (G) 


2 2 

^ 00} ) y <|>Y  y 

A-»  A-,  A *) 

i Z 12 


(Ow)y  • = 4 y - , 

h"2  -1a2  1200  1 

and  for  (9) , approximately 


C-. .'i 2 a ,2  13  rn  . a2,2  f16  ,2.  2 2 

7 ^XXX2  ' 1200  [fl  13}  (35  ct2)  ]a  al 


ly  \35  2J 


L. 


0 


FARLIE,  D.  J.  G.  (1960).  The  performance  of  some  correlation  coefficients 
for  a general  bivariate  distribution,  Biometrika , 47,  307-323. 

JOHNSON,  N.L.  ar.d  X0T7. , 7.(1976).  On  some  generalized  Farlie-Gumbel-I brgenstem 
distributions  — II.  Regression,  correlation  and  further  generalizations. 
Mineo  Series  No.  1030,  Institute  of  Statistics,  University  of  North 
Carolina.  (A  slightly  different  version  will  appear  in  Communications 
in  Statistics.) 

A 

REI-JYI , A.  (1970).  Probability  Theory.  Anster lam : North  Holland  Publishing  Co. 

SCHi/IilZER,  B.  u NOLFF,  E.  F.  (1970).  Sur  nne  insure  de  d^oenaance  pour  les 
variables  aleatoires,  C.P:  Acad.  Soi. , Paris , Sir.  A,  235,  659-661. 


UNCLASSIFIED 


EC.URITV  CLASSIFICATION  OF  THIS  PAGE  rWhu.  Pmlm  Enfrtd) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


T REPORT  NUMBER 


2 GOVT  ACCESSION  NO. 


PIENT'S  CATALOG  NUMBER 


* . T4TLE  (and  SuhtttU) 


/Tlependence 


bivariate 


'properties  of  iterated  generalized 
Farl ie-Gumbel -Morgenstern  distributions^ 


E OF  REPORT  &.PEJUOO. COVERED 

JECHNICAL  / , j 

^umbeIr"" 


6.  PERFORMING  ORG.- -REPORT  fjujjj 

Mimeo  Series  No  .If  0^ 


■p.  aothor?.)-  w . 

^Sy'Kotz  St  N.  L./Gohnson  j 


».  contract  or 


VS) 

1 


l)AA  G29-74-C-0Q3O 


MBERf*) 


y. 


IQ.-  PBOCnAM  BLrMPNT-FTOJECT~.~Tr?R' 


S.  PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 

Department  of  Statistics 
University  of  North  Carolina 
Chapel  Hill,  North  Carolina  27514 


AREA  4 WORK  UNIT  NUMBERS 


It.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

U.  S.  Army  Research  Office 
Box  1221 

Research  Triangle  Park,  NC  27709 


7 7} 


12.  HEPSrt  daYF 

Mar4fc  077 


'WJMBgR  OF  P AO 

10 


l 


M.  MONITORING  AGENCY  NAME  5 AODRESS/T/  dlllerent  from  Controlling  Oltlce) 

i / AV.  ’f./j’  1 v | 

~t_r-  4 ' 1 


15.  SECURITY  CLASS,  (cl  thle  report) 

UNCLASSIFIED 


15«.  DECL  ASSI  FI  CAT!  On/  DOWNGRADING 
SCHEDULE 


16-  DISTRIBUTION  STATEMENT  (of  thle  Report) 

Approved  for  Public  Release:  Distribution  Unlimited 

( TsJJZp,  i 


17.  DISTRIBUTION  ST  AT  EM  pH  T -fir 


Hte,  ft  iffffilWH  IfOOl  Report) 


18.  SUPPLEMENTARY  NOTES 


*19  KEY  WORDS  (Continue  on  reveree  aide  II  neceeemry  end  Identify  by  block  number) 

j 

f 

! Measures  of  Dependence,  Bivariate  Distributions,  Generalized  Farl ie-Gumbel - 
Morgenstern  Distributions 


20  ABSTRACT  (Continue  on  reveree  aide  It  neceeemry  end  Identity  by  block  number) 

Farl ie-Gumbel -Morgenstern  distributions  allow  for  only  a limited  degree  of 
dependence  between  the  variables.  Some  modifications,  designed  to  increase 
dependence,  are  discussed,  ft  is  found  that  only  a limited  increase  in 
dependence  is  feasible  with  these  modifications. 


DD  , 1473  EOITION  OF  I NOV  ••  IS  OBSOLETE 

UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PACE  (Whan  Data  Bnltmd) 


